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Abstract 

The  main  objective  of  this  research  project  is  to  demonstrate  through  a  combination  of  experimental 
measurements  and  numerical  simulations  that  ultrahydrophobic  surfaces  can  be  used  to  delay  the 
transition  to  turbulence  and  dramatically  reduce  drag  in  both  external  and  internal  turbulent  flows. 
These  enhancements  could  have  a  profound  effect  on  a  huge  variety  of  existing  technologies,  resulting 
in  benefits  ranging  from  a  reduction  in  the  pressure  drop  in  pipe  flows,  to  improved  microfluidic  devices, 
to  an  increase  in  range  of  small  submersible  vehicles.  In  short,  almost  any  device  that  uses  or  interacts 
with  liquids  has  the  potential  to  benefit  from  the  ultrahydrophobic  surface  technologies  that  will  be 
developed  through  the  research  performed  for  this  contract. 

Technical  Section 

See  pages  4-16 

Final  Progress  Statement 

Over  the  past  year  we  have  been  investigating  the  turbulent  drag  reduction  properties  of 
ultrahydrophobic  surfaces  through  both  detailed  experiments  and  direct  numerical  simulations.  Our 
results  demonstrate  drag  reduction  of  up  to  75%  using  ultrahydrophobic  surfaces  with  slip  lengths  of 
over  120  microns.  The  drag  reduction  is  observed  to  begin  after  the  flow  has  become  fully  turbulent  at 
or  near  Re  =  4000  and  to  increase  with  increasing  Reynolds  number  and  microfeature  spacings.  The 
simulation  results  were  recently  published  in  the  Journal  of  Fluid  Mechanics  with  a  second  paper 
submitted  to  Physics  of  Fluids.  Our  experimental  results  were  recently  published  in  Physics  of  Fluids  and 
my  recent  Annual  Review  of  Fluid  Mechanics  article. 
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In  May  of  2009  we  applied  for  a  provisional  patent  based  on  our  ONR  funded  research.  The  patent  is 
entitled  "Design  of  Superhydrophobic  Surfaces  for  Drag  Reduction."  We  hope  to  license  and 
commercialize  this  technology  in  the  coming  year. 
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Technical  Objectives 

The  interaction  between  solid  surfaces  and  liquids  is  of  fundamental  importance  in  engineering  flows 
where  solid  surfaces  are  the  primary  means  for  controlling  or  manipulating  fluids.  We  will  show  that  by 
treating  a  solid  surface  to  make  it  ultrahydrophobic  it  is  possible  to  significantly  enhance  how  the  surface 
interacts  with  a  flowing  liquid.  In  particular,  in  our  recent  work,  we  have  shown  that  ultrahydrophobic 
surfaces  can  be  utilized  to  reduce  drag  [2,  3]  and  enhance  mixing  [4]  in  low  to  moderate  Reynolds 
number  internal  flows.  The  main  objective  of  this  research  project  is  to  demonstrate  through  a 
combination  of  experimental  measurements  and  numerical  simulations  that  ultrahydrophobic  surfaces 
can  also  be  used  to  delay  the  transition  to  turbulence  and  dramatically  reduce  drag  in  both  external  and 
internal  turbulent  flows.  These  enhancements  could  have  a  profound  effect  on  a  huge  variety  of  existing 
technologies,  resulting  in  benefits  ranging  from  a  reduction  in  the  pressure  drop  in  pipe  flows,  to 
improved  microfluidic  devices,  to  an  increase  in  range  of  small  submersible  vehicles.  In  short,  almost 
any  device  that  uses  or  interacts  with  liquids  has  the  potential  to  benefit  from  the  ultrahydrophobic 
surface  technologies  that  will  be  developed  through  the  research  proposed  here. 

The  presence  of  drag  reduction  using  ultrahydrophobic  surfaces  had  been  suggested  previously  by  a 
number  of  other  authors  [5-7].  However,  prior  to  our  experiments  and  numerical  simulations  [2-4],  the 
results  of  other  researchers  were  somewhat  ambiguous  and  a  number  of  different  theories  were  presented 
to  explain  the  observed  drag  reduction.  In  our  previous  work,  a  detailed  experimental  study  was 
performed  which  unequivocally  demonstrated  the  existence  of  drag  reduction  in  laminar  flows  with 
ultrahydrophobic  surfaces.  Ultrahydrophobic  surfaces  are  engineered  by  taking  materials  with  micron  or 
nanoscale  surfaces  roughness  and  chemically 
treating  them  to  make  them  hydrophobic.  When 
water  is  brought  in  contact  with  an 
ultrahydrophobic  surface,  the  water  is  completely 
supported  by  the  peaks  in  the  surface  topology  as 
seen  in  Figure  1.  Because  of  the  hydrophobicity 
of  these  microscale  and  nanoscale  protrusions, 
the  water  does  not  move  into  the  pores  on  the 
surface,  rather  it  remains  in  contact  with  only  the 
peaks  of  the  surface  topology  resulting  in  a 
shear-free  air-water  interface.  The  underlying 
physical  mechanism  for  drag  reduction  is  a  slip 
along  the  shear-free  air-water  interface  supported  Figure  1.  Micrographs  of  a  water  droplet  in 

between  the  peaks  of  micro  or  nanoscale  contact  with  an  ultrahydrophobic  lotus  leaf.  The 

protrusions  present  on  the  ultrahydrophobic  droplet  is  supported  by  the  hydrophobic  micro-posts 

populating  the  surface  of  the  lotus  leaf  (from  [1]) . 
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surface.  This  physical  mechanism  was  validated  for  laminar  flows  through  direct  observation  of  the 
interface  using  laser  confocal  metrology,  direct  velocity  measurement  using  micro  particle  image 
velocimetry  (p-PIV),  and  direct  comparison  to  numerical  simulations  and  analytical  theory  [2-4].  One  of 
the  major  objectives  of  this  project  is  to  demonstrate  that  these  same  surfaces  can  also  reduce  drag  when 
the  flow  is  turbulent. 

Once  we  have  fully  understanding  the  physical  mechanism  of  ultrahydrophobic  drag  reduction  our 
goal  will  be  to  configure  the  microscopic  features  of  the  surface  to  optimize  its  performance  under 
different  design  criteria.  We  originally  used  etched  silicon  wafers  in  order  to  precisely  control  the  micro¬ 
geometry  and  therefore  unambiguously  control  the  micro-characteristics  of  and  flow  characteristics  near 
the  ultrahydrophobic  surface,  but  we  have  recently  been  able  extend  our  technique  to  more  flexible 
materials  like  polydimethylsiloxane  (PDMS)  through  soft  lithography  and  pressed  Teflon.  The  use  of 
these  well-designed  surfaces  has  and  will  enable  us  to  perform  exactly  matching  computational  fluid 
dynamics  of  the  flow  near  these  ultrahydrophobic  surfaces  [3,  4], 

In  summary,  the  technical  objectives  of  this  research  project  are  to:  (1)  demonstrate  that  turbulent 
drag  reduction  is  possible  with  these  surfaces;  (2)  confirm  that  the  transition  to  turbulence  is  delayed;  and 
finally  to  (3)  develop  techniques  for  fabricating  large  area  sheets  of  ultrahydrophobic  surfaces  with 
nanometer  length  features  for  use  in  deep-sea  applications.  The  goal  is  to  understand  precisely  how  these 
materials  behave,  maximize  their  efficiency  and  determine  under  what  circumstances  they  are  the  most 
useful. 


Technical  Approach 

Experiments 

An  experimental  flow  cell  was  fabricated  for  taking  pressure  drop  and  particle  image  velocimetry 
(PIV)  measurement  of  the  turbulent  flow  past  ultrahydrophobic  surfaces.  An  image  of  the  flow  cell  is 
shown  in  Figure  2.  In  the  flow  cell  the  top  wall  is  always  smooth  and  optically  clear  for  PIV  and  the 
bottom  wall  is  changeable  and  could  be  smooth  or  ultrahydrophobic.  The  flow  cell  is  2in  wide,  3ft  long 
and  has  a  height  that  is  variable  from  l/8in  to  l/2in  high.  A  gravity  feed  system  is  used  to  isolate  the 
flow  from  the  vibrations  and  fluctuations  of  the  pump  which  is  capable  of  driving  flows  up  to  a  Reynolds 
number  of  approximately  Re  =  UH/v  =  10,000.  Here  U  is  the  average  velocity  in  the  channel,  H  is  the 
height  of  the  channel  and  v  is  the  viscosity  of  the  fluid.  A  flow  meter  is  used  to  determine  the  flow  rate 
and  average  Reynolds  number  in  the  channel  and  a  series  of  differential  pressure  transducers  are  used  to 
measure  the  pressure  drop.  Additionally,  a  PIV  system  is  installed  using  an  Argon  Ion  laser  and  a  high¬ 
speed  video  camera  to  take  instantaneous  and  time  averaged  measurements  of  the  velocity  and  turbulent 
statistics  within  the  flow.  The  images  are  processed  using  the  commercial  PIV  software  LaVision. 


Figure  2,  Experimental  flow  cell  used  for  ultrahydrophobic  drag  reduction  measurements. 


5 


Ultrahydrophobic  surfaces  created  by  using 
either  standard  photolithographic  techniques  to  etch 
silicon  wafers  [2,  8]  or  soft  lithography  to  produce 
precise  micropattems  in  PDMS  [9]  or  through  a 
new  technique  we  have  developed  to  cold-work 
metal  mesh  screens  into  Teflon.  A  schematic 
diagram  outlining  the  fabrication  procedure  for 
lithographically  etched  ultrahydrophobic  surfaces 
is  shown  in  Figure  3.  First,  a  CAD  package  is  used 
to  design  the  size,  spacing  and  alignment  of  the 
desired  micro-patterned  surface.  The  file  is  then 
printed  on  a  high-resolution  transparency  at  5080 
dpi.  The  transparency  serves  as  a  mask  for  contact 
photolithography  using  a  positive  photoresist  on  a 
silicon  wafer  [10].  At  this  point  there  are  two 
options;  to  use  either  wet  or  DRIE  etching 
techniques  to  create  the  desired  pattern  in  the  silcon 
surface  as  shown  in  Figure  3  or  to  use  a  photoresist 
like  SU-8  which  can  create  robust  features  without 
the  need  for  etching  that  can  be  up  to  100pm  tall. 

The  silicon  wafer  is  in  turn  used  as  a  master  for 
creating  a  negative  of  its  pattern  in  PDMS  using 
soft  lithography  [9].  The  PDMS  surface  is  then 
cleaned  using  a  plasma  etcher  and  silanized  under  a 
vapor  reaction  with  an  organosilane  reagent  [11]. 

These  ultrahydrophobic  surfaces  can  have 
advancing  contact  angles  as  large  as  ^=165°.  Soft 
lithography  is  non-photolithographic  technique  developed  for  rapidly  replicating  a  micropattem.  The 
micropattems  of  interest  are  transferred  to  the  PDMS  through  replica  molding  from  a  master  creating  a 
negative  replica  in  the  PDMS.  The  master  can  be  a  photolithographically  etched  silicon  wafer  or  any 
other  surface  with  the  desired  surface  roughness.  After  being  spun  onto  or  compressed  into  the  surface 
of  the  master,  the  PDMS  is  cured  in  an  oven  and  the  replica  is  peeled  from  the  master  creating  a 
micropatterned  surface  on  a  flexible  elastomer  substrate.  It  is  naturally  hydrophobic  without  further 

chemical  reaction  and  it  can  be  applied  to 
almost  any  substrate  or  scaffold  while 
maintaining  the  ordered  and  precise  patterning 
that  makes  the  use  of  silicon  wafers  so 
desirable. 

To  fabricate  a  very  large 
ultrahydrophobic  surface,  a  number  of  6in 
PDMS  surfaces  were  patched  together  using 
fast  curing  PDMS  as  glue.  Great  care  was 
taken  to  minimize  the  profile  of  the  seams.  The 
quality  of  the  patches  were  confirmed  with 
surface  profilometry  measurements.  We  are 
also  working  on  a  number  of  techniques  to 
produce  long  seemless  ultrahydrophobic 
surfaces.  We  are  currently  developing  a  roll- 
to-roll  system  capable  of  transferring  precise 


Figure  4.  Example  of  an  ultrahydrophobic  surface 
made  from  PDMS  with  30|im  wide  microridges 
spaced  30|im  apart 


Silicon  dioxide 
Silicon 


Photolithography 


Silicon  dioxide 
etching 


Silicon  etching 
Silanized  area 


Silanization 


Figure  3,  Schematic  diagram  of  the  method 
used  for  creation  of  ultrahydrophobic  surfaces. 
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Figure  5.  Schematic  diagram  of  an  ultrahydrohobic  microridge  currently  being  modeled.  Note  that 
the  ridge  height  does  not  exist  in  code.  The  ‘top’  of  the  ridges  have  a  no-slip  boundary  condition, 
and  the  ‘space’  between  the  ridges  have  a  slip  boundary  condition. 


patterns  onto  a  film  in  a  continuous  and  seamless  fashion. 

Experiments  have  been  performed  with  a  range  of  different  surfaces.  For  the  purpose  of  this 
progress  report  we  will  focus  on  the  results  from  two  surfaces:  a  30)im  wide  microridges  spaced  30|xm 
apart  as  seen  in  Figure  4  and  a  60)im  wide  microridges  spaced  60)im  apart  for  which  we  have  the  most 
data,  but  will  also  discuss  surfaces  with  other  geometries  as  well. 

Numerical  Simulations 

Direct  numerical  simulations  are  currently  being  performed  for  flow  over  a  series  of  model 
ultrahydrophobic  surfaces  chosen  to  directly  match  the  surfaces  used  in  the  experiments.  In  the 
simulations,  the  top  of  the  microposts  or  microridges  are  assumed  to  be  no-slip  and  the  air-water 
interface  between  them  is  assumed  to  be  shear-free  and  flat.  Both  of  these  assumptions  are  reasonable, 
however,  under  certain  conditions,  the  air-water  interface  can  deflect  and  the  recirculation  of  air  within 
the  gaps  between  the  microposts  or  microridges  can  produce  some  drag  along  the  air-water  interface. 
Quantifying  the  contribution  of  both  of  these  effects  is  important  and  is  planned  for  the  final  year  of  this 
project.  The  simulation  domain  currently  contains  ten  or  more  microridges  and  periodic  boundary 
conditions  on  the  side  walls  as  seen  in  Figure  5.  The  numerical  algorithm  used  to  simulate  the  has  been 
under  development  for  almost  a  decade  [12,  13]  and  over  the  last  year  it  has  been  modified  to  apply 
directly  to  turbulent  flows  past  ultrahydrophobic  surfaces.  The  code  uses  non-uniform  grid  spacing  to 
place  more  resolution  near  the  channel  walls  and  a  Cartesian  staggered  mesh  method  [14]  with  exact 
projection  [15]  for  the  pressure  solution.  The  Navier-Stokes  equations  are  solved  using  a  conjugate 
gradient  solver  and  a  modified  Runge-Kutta  time  marching  method  [16].  The  code  is  fully  parallelized 
[17]  (with  MPI)  for  PC  cluster  architectures.  The  code  is  currently  running  on  the  Navy’s  supercomputer 
sight  in  Alaska  with  resolutions  of  128^,  256^  and  512^. 


Progress  Statement  Summary 

Over  the  course  of  the  past  year  we  have  been  investigating  the  turbulent  drag  reduction  properties  of 
ultrahydrophobic  surfaces  through  both  detailed  experiments  and  direct  numerical  simulations.  On  the 
experimental  side,  we  have  fabricated  a  flow  cell  and  developed  a  number  of  techniques  for  producing 
large-scale  surfaces  containing  precise  micron-scale  patterns  of  microridges  and  microposts.  Using  both 
pressure  drop  measurements  and  particle  image  velocimetry  measurements  our  results  demonstrate  drag 
reduction  of  up  to  75%  using  ultrahydrophobic  surfaces  with  slip  lengths  of  over  120)im.  The  drag 
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reduction  is  observed  to  begin  after  the  flow  has  become  fully  turbulent  at  or  near  Re  =  4000  and  to 
increase  with  increasing  Reynolds  number  to  the  limit  of  our  experiments  at  Re  =  10,000.  On  the 
numerical  simulation  side,  we  modified  an  existing  direct  numerical  simulation  code  to  perform 
simulations  with  walls  having  periodic  slip  no-slip  boundary  conditions  to  simulate  the  flow  over 
hydrophobic  microridges  and  microposts.  The  code  is  now  fully  operational  and  parallelized  having 
been  benchmarked  against  simulation  results  in  the  open  literature  for  flow  within  rectangular  channel 
with  smooth  walls  at  128^,  256^  and  512^.  We  have  simulating  turbulent  flow  over  ultrahydrophobic 
surfaces  containing  both  microridges  and  microposts  that  match  and  extend  those  used  in  the 
experiments.  The  simulations  agree  quite  well  with  the  experimental  results  and  demonstrate  drag 
reductions  and  slip  velocities  which  increase  with  increasing  spacing  between  microfeatures  and 
Reynolds  number.  The  DNS  work  for  128^  simulations  was  recently  published  in  the  Journal  of  Fluid 
Mechanics  [18]  as  a  fast  track  article  and  a  follow  up  article  with  results  for  higher  Reynolds  numbers 
was  recently  submitted  to  Physics  of  Fluids  [19].  Our  experimental  results  were  recently  accepted  as  a 
featured  article  in  Physics  of  Fluids  [20]  and  a  second  paper  with  the  results  from  a  number  of  different 
surface  geometries  is  currently  being  prepared  for  submission  to  Journal  of  Fluid  Mechanics  [21].  This 
work  has  been  presented  at  a  number  of  national  and  international  conferences  and  workshops  and  has 
been  featured  in  my  recent  Annual  Review  of  Fluid  Mechanics  article  [22]  entitled  “Slip  on 
Superhydrophobic  Surfaces.” 


Progress 

Experiments 

To  date  we  have  concentrated  our  efforts  on  taking  measurements  of  the  flow  past  a  smooth 
PDMS  surface  and  a  series  of  ultrahydrophobic  surfaces  consisting  of  lOpm,  15p,m,  30p,m  and  60p,m 
wide  microridges  spaced  between  15p,m  and  60p,m  apart.  For  sake  of  brevity,  in  this  progress  report,  we 
will  focus  primarily  on  two  geometries  a  30p,m  wide  microridges  spaced  30p,m  apart  as  seen  in  Figure  4 
and  a  60pm  wide  microridges  spaced  60pm  apart,  but  we  will  also  present  highlights  from  the  best 
surfaces  studied  to  date. 

A  typical  set  of  velocity  profiles,  resulting  from  PIV  near  the  superhydrophobic  wall  for  the  60-60 
ridge  surface  is  shown  in  Figure  6a  for  a  range  of  Reynolds  number  between  2700<  Re  <  8200.  The 
effect  of  the  superhydrophobic  wall  is  not  observed  for  the  low  Reynolds  number  experiments.  At  the 
low  Reynolds  numbers,  the  turbulent  velocity  profiles  just  past  transition  are,  to  the  limit  of  our 
measurements,  equivalent  to  smooth  profiles  at  identical  Reynolds  numbers  This  is  not  unexpected  for 
the  data  points  in  the  laminar  or  transitional  regime  [2,  3].  In  laminar  flows,  superhydrophobic  surfaces 
of  similar  size  and  geometry  demonstrated  slip  lengths  which  were  independent  of  flow  rate  and 
approximately  b  =  25pm  [3].  For  pressure  driven  flow  between  two  infinite  parallel  plates  separated  by  a 
distance  2S  the  volume  flow  rate  per  unit  depth  is  given  by 


26^ 

fl  b  I 

1  ...  _ 

_3  b  +  2S^ 

For  a  given  pressure  gradient,  dpidx,  and  fluid  viscosity,  p,  the  volume  flow  rate  can  be  significantly 
enhanced  only  if  the  slip  length  is  comparable  to  the  channel  height.  Previous  laminar  regime  studies 
over  similar  superhydrophobic  microfeatures  measured  slip  lengths  of  =  25pm  independent  of 
Reynolds  number  [2].  In  our  channel  geometry,  such  laminar  flow  slip  lengths  would  produce  a  drag 
reduction  of  around  1%.  Additionally,  for  small  slip  lengths,  the  expected  slip  velocity  can  be 
approximated  by  =AUbl 5  which  should  also  be  on  the  order  of  only  a  couple  of  percent  of  the 

average  free  stream  velocity,  U,  and  below  the  resolution  of  our  PIV  measurements.  As  the  Reynolds 
number  is  increased  and  the  flow  becomes  fully  turbulent,  however,  a  substantial  slip  velocity,  and  slip 
lengths  greater  than  b  -  25pm,  are  observed  along  the  superhydrophobic  wall.  The  presence  of  an  air 


8 


Figure  6.  a.)  Velocity  profiles  over  a  microridge  surface  w  =  60pm  d  =  60pm  showing  the  development 
of  significant  slip  velocities  with  increasing  Reynolds  number  from  2700  (A)  to  8200  (■).  (Inset) 
Velocity  profiles  near  the  wall  demonstrating  prominent  slip  velocities.  Reynolds  numbers  are:  2700 
(A),  3900  (A),  4840  (0),  5150  (♦),  6960  (□),  8200  (■).  For  clarity,  the  modified  Spalding  fits  (— )  from 
Equation  3  are  only  overlaid  on  the  profiles  corresponding  to  /?e=2700  and  /?e=8200.  b.)  Velocity 
profiles  over  the  w  =  30pm  d  =  30pm  microridge  surface  demonstrate  slip  velocity  behavior  consistent 
with  that  observed  on  the  60-60  surface,  but  reduced  in  magnitude.  Reynolds  numbers  range  from  4970 
(o)  to  7930  (V).  Larger  feature  spacing  performs  better  for  a  given  Reynolds  number.  Reynolds 
numbers  are:  4970  (o),  5400  (♦),  6800  (A),  7160  (■),  7930  (V)  The  modified  Spalding  fits  (— )  are 
overlaid  on  the  profile  corresponding  to  /?e=7930. 
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water  interface  is  visually  apparent  on  the  superhydrophobic  surface  giving  it  a  silvery  appearance.  This 
result,  due  to  the  differing  indices  of  refraction  and  slight  curvature  of  the  interface,  was  observed 
throughout  the  range  of  testing  giving  us  confidence  that  the  interface  was  maintained  for  all  of  the 
experiments  reported  in  this  paper. 

As  the  inset  of  Figure  6a  clearly  shows,  the  magnitude  of  the  slip  velocity  was  found  to  increase 
with  increasing  Reynolds  number.  Similar,  although  less  pronounced,  trends  were  observed  for  the  30-30 
ridge  case  as  seen  in  Figure  6b.  Significant  deviation  from  no-slip  behavior  is  noted  past  a  Reynolds 
number  of  approximately  Re  =  4000  for  both  the  30-30  and  60-60  ridged  cases.  Above  these  Reynolds 
numbers,  a  nearly  linear  increase  in  the  slip  velocity  with  increasing  Reynolds  number  was  observed  for 
each  of  the  superhydrophobic  surfaces  used.  A  maximum  slip  velocity  of  nearly  40%  the  mean  channel 
velocity,  lU  =  0.4  was  observed  for  the  60-60  ridged  case  at  the  highest  Reynolds  numbers  tested. 


In  order  to  determine  both  the  shear  stress  and  slip  velocity  at  the  smooth  and  superhydrophobic 
walls,  the  PfV  velocity  fields  were  fit  to  a  modified  Spalding  equation  for  turbulent  velocity  profile 
above  a  flat  plate  [23], 


~<»p) 


+  e 


-205 


-0.41 


The  Spalding  equation  is  an  empirical  fit  to  experimental  turbulent  velocity  profile  data  that  covers  the 
entire  wall  region  through  the  log  layer  [24].  This  allows  the  fit  to  be  applied  farther  into  the  channel,  to 
determine  the  wall  shear  stress  more  accurately  using  a  greater  number  of  data  points  than  would  be 
available  within  the  viscous  sublayer.  Wall  shear  stress  enters  the  equation  in  the  definition  of  the 
velocity,  and  position  y^,  in  wall  units.  To  account  for  slip,  each  instance  of  the  velocity  in  wall 
units,  =Uyjp/r^  ,  in  the  Spalding  equation  was  replaced  by  the  difference  The  fit  was 


performed  by  a  numerical  routine  given  an  initial  value  for  slip  velocity  extrapolated  from  a  coarse  linear 
fit  of  near  wall  data  points.  An  initial  wall  shear  stress  was  determined  by  minimizing  the  error  in  the  fit. 
Subsequent  iterations  were  performed  on  wall  slip  velocity  and  wall  shear  stress  to  minimize  the  standard 
error  of  the  fit  over  the  interval  0  <  50.  The  resulting  fits  were  accurate  to  better  than  4%  at  a  95% 

confidence  interval.  The  results  were  not  appreciably  different  if  the  fit  is  taken  to  y^  =  1 00.  The  size  of 
the  PIV  correlation  window  was  chosen  to  be  0.2mm,  For  the  frame  rates  used,  the  resulting  particle 
displacements  within  the  correlation  window  were  typically  much  less  than  25%  of  the  window  in  the 
viscous  sublayer  and  less  than  33%  of  the  window  everywhere  for  Reynolds  numbers  less  than  Re  < 
4500.  Large  particle  displacements  were  observed  far  from  the  wall  at  the  highest  Reynolds  numbers, 
however,  no  noticeable  effects  were  observed  on  the  resulting  profiles. 

As  seen  in  Figure  6,  the  resulting  fits  of  Equation  3  to  the  velocity  profiles  are  excellent  with  and 
without  slip,  which  instills  confidence  in  the  values  of  shear  stress  calculated  from  the  velocity  gradient 
extrapolated  to  the  wall,  jj  .  The  maximum  slip  velocity  and  observed  wall  shear  stress 


reductions  correspond  to  slip  lengths  of  b>  TOpm  for  the  30-30  microridges  and  b  >  120pm  for  the  60-60 
microridges.  Larger  slip  velocities  and  slip  lengths  were  measured  for  turbulent  flow  past 
superhydrophobic  surfaces  with  larger  microfeature  spacings  even  as  the  percentage  of  shear-free 
interface  was  kept  constant  at  w/(w^d)  -  0.5,  as  has  been  observed  in  the  laminar  flow  measurements 
over  superhydrophobic  surfaces  [3].  This  observation  is  consistent  previous  laminar  flow  studies[2,  3] 
and  with  the  predictions  of  DNS  in  turbulent  flows  [25].  Additionally,  Ybert  et  ai  [26]  showed  through 
a  scaling  argument  that  in  laminar  flows  one  expects  the  slip  length  to  scale  linearly  with  the 
microfeature  spacing,  b<x:(w  +  ci). 

In  Figure  7,  direct  measurements  of  the  pressure  drop  per  unit  length  of  channel,  dpH,  are  shown 
for  a  smooth  PDMS  surface  and  the  superhydrophobic  surface  containing  60pm  ridges  spaced  60pm 
apart  and  a  surface  containing  10pm  ridges  spaced  60pm  apart  in  an  identical  channel.  The  result 
predicted  by  the  Colebrook  equation  [27]  for  a  perfectly  smooth  channel  of  the  same  dimension  is  plotted 
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for  reference.  The  pressure  drop  per  unit 
length  is  directly  related  to  the  channel 
geometry  and  the  wall  shear  stress, 
dpH  -t^{\^25IW^I  5  ^  so  it  provides  a 

second  method  for  measuring  drag 
reduction.  Significant  drag  reduction  is 
initially  noted  by  a  leveling  off  of  the  in 
the  pressure  drop  during  the  transition  from 
laminar  to  turbulent  flow  between 
Reynolds  numbers  of  2000  <  Re  <  3000. 

These  data  indicate  a  delay  in  the  transition 
to  fully-developed  turbulent  flow  and  an 
increased  drag  reduction  with  increased 
microridge  spacing.  Additionally,  for 
Reynolds  numbers  greater  than  Re  >  3000 
the  pressure  drop  over  the  60-60 
surperhydrophobic  surface  grows  at 
roughly  half  the  rate  of  pressure  drop  over 
the  smooth  surface  while  the  60-10  grows 
at  less  than  one  quarter  the  rate  of  the 
smooth  surface.  The  Colebrook  line, 
accurately  fits  the  turbulent  flow  data  from 
the  smooth  surface,  and  the  predicted 
laminar  flow  result  passes  through  the  microridge  data  in  the  laminar  region  below  Re  <  2200.  This 
result  is  consistent  with  those  predicted  by  Equation  2  and  observed  by  PIV.  As  noted  before  there  is  no 
measurable  drag  reduction  or  slip  velocity  for  the  present  channel  geometry  in  the  laminar  regime. 

In  Figure  8,  the  wall  shear  stresses,  calculated  from  the  Spalding  fit  to  the  velocity  profiles 
and  from  pressure  measurements  of  smooth,  30-30  and  60-60  channels  are  non-dimensionalized  to  form  a 
coefficient  of  friction,  Cj^  =  2r^^/,  /  pU^ ,  and  plotted  as  a  function  of  Reynolds  number.  For  comparison, 

the  Colebrook  prediction  of  friction  coefficient  for  the  present  perfectly  smooth  channel  is  superimposed 
over  the  data  in  Figure  8.  Friction  coefficient  was  selected  to  account  for  small  variations  in  channel 
height  existing  between  the  pressure  drop  and  PIV  experiments.  As  previously  indicated,  the  friction 
coefficients  of  the  smooth  wall,  calculated  from  PIV,  and  that  of  the  smooth  channel,  determined  from 
pressure  drop,  are  in  good  agreement  with  each  other  as  well  as  with  the  Colebrook  prediction.  At  low 
Reynolds  numbers,  in  the  absence  of  any  quantifiable  slip  at  the  superhydrophobic  wall,  the  coefficient 
of  friction  for  all  cases  tracks  with  that  of  the  smooth-walled  channel.  At  larger  Reynolds  numbers, 
where  slip  velocities  are  observed,  the  coefficients  of  friction  of  the  superhydrophobic  surfaces  were 
found  to  lie  well  below  those  of  the  smooth  channels.  The  drag  reduction  was  found  to  increase  with 
increasing  Reynolds  number  and  feature  spacing,  becoming  more  significant  for  Re  >  5000  as  observed 
in  the  pressure  measurements.  The  60-10  surfaces  appear  to  fully  relaminarize  the  flow  resulting  in 
nearly  75%  drag  reduction  at  high  Reynolds  numbers.  The  PIV  measurements  of  the  channel  with  a  30- 
30  superhydrophobic  microridge  surface  on  one  wall  and  a  smooth  no-slip  surface  on  the  opposing  wall 
show  a  somewhat  smaller  drag  reduction  than  that  which  is  noted  by  pressure  drop  along  with  two 
superhydrophobic  walls.  This  result  is  likely  due  to  differences  in  the  flowcell  geometry,  specifically, 
the  presence  of  the  smooth  wall  in  the  PfV  measurements,  which  was  necessary  to  have  transparency  for 
flow  visualization.  The  smooth  wall  has  a  higher  wall  shear  stress  than  the  superhydrophobic  surface 
resulting  in  an  asymmetric  velocity  profile  and  an  increase  in  the  turbulence  intensity  near  the  smooth 
wall.  These  observations  were  also  made  by  Martell  et  al.  [25]  [28]  for  a  DNS  of  channel  flow  with  a 
single  superhydrophobic  wall.  Observed  drag  reductions  and  slip  velocities  are  in  good  agreement  with 


Figure  7.  Pressure  drop  measurements  for  flow  through  a 
rectangular  channel  with  a  smooth  walls  (A)  and  with  two 
walls  containing  superhydrophobic  microridges  with 
>v=60pm  and  ^/=60pm  (■)  and  >v=60pm  and  c/=10pm  (0). 
The  Colebrook  line  ( — )  is  shown  for  a  smooth  channel. 
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Figure  8.  Coefficient  of  friction  for  various  surfaces  calculated 
from  both  PIV  and  pressure  measurements.  Smooth  surfaces 
(A)  and  superhydrophobic  surfaces  containing  w=30//m  wide 
microridge  spaced  c/=30ywm  apart  (•)  are  shown  for  PIV 
measurements  of  a  channel  with  a  single  superhydrophobic 
wall.  Pressure  drop  measurements  from  channels  with  two 
smooth  walls  and  two  superhydrophobic  walls  containing 
>v=30p.m  and  (#=30pm  microridges  (o),  w=60|im  ^i=60|im 
microridges  (■)  and  >v=60pm  rf=10pm  microridges  (*)  are  also 
shown.  The  predictions  of  the  friction  coefficient  for  a  smooth 
channel  is  also  shown  ( — )  in  both  the  laminar  and  turbulent 
regimes.  Transition  occurs  around  Re  =  2100. 


predictions  for  a  DNS  at  Re^  =180, 

corresponding  to  an  experimental 
/?e=5300  in  the  PIV  data.  DNS 
slightly  over  predicts  slip  velocity, 
and  slightly  under  predicts  drag 
reduction  at  1 1%  and  reports 
enhanced  performance  with 
increasing  microfeature  size,  as 
observed  in  the  experiments.  It 
should  also  be  noted  that  DNS  of 
Martell  et  al  [25]  [28]  does  not 
include  interface  deflection  or 
compliance  effects.  Drag  reduction 
calculated  from  PIV  data  are  in 
excellent  agreement  with  the  slip 
length  boundary  condition  DNS  of 
Min  and  Kim  [29]  and  predictions  of 
Fukagata  et  al.  [30]  for  stream  wise 
slip.  Both  groups  reported 
approximately  21%  drag  reduction 
[29,  30]  at  the  same  dimensionless 
slip  length  and  friction  Reynolds 
number  observed  in  the  present 
experiments  at  /?e=5300.  Given  the 
challenges  of  directly  matching  DNS 
and  experiments,  these  results  are 
quite  encouraging. 

The  turbulent  drag  reduction, 

Dh  =  {^^no-sllp  -^Sh)/  ^no-,Hp  >  ^38 

computed  as  the  percent  difference 
in  shear  stress  at  the 
superhydrophobic  and  no-slip  wall 


and  is  presented  in  Figure  9  as  a 
function  of  Reynolds  number.  Drag  reduction  is  presented  rather  than  slip  length  because  the  slip  length 
is  difficult  to  quantify  from  the  pressure  drop  measurements  in  turbulent  flows.  The  slip  length 
calculated  from  PIV  data  is  insignificant  in  the  laminar  region  and  obtains  a  maximum  value  greater  than 
b  =  70pm  for  30-30  and  greater  than  b  =  120pm  for  60-60  ridges.  In  the  present  experiments,  a  maximum 
drag  reduction  of  approximately  50%  was  observed  for  both  microridge  geometries  once  a  suitably  high 
Reynolds  number  was  achieved.  Drag  reduction  is  initiated  at  a  critical  Reynolds  number  in  the  turbulent 
regime.  For  the  microridges  under  present  consideration,  the  critical  Reynolds  number  was  determined 
to  be  Recrit  ~  2500.  This  Reynolds  number  is  at  or  just  past  the  transition  to  turbulent  flow.  This 
observation,  along  with  the  noted  lack  of  drag  reduction  in  the  laminar  regime,  suggest  that  the 
underlying  physical  cause  of  the  observed  turbulent  drag  reduction  must  relate  to  the  unique  structure  of 
wall-bounded  turbulent  flow. 

The  physical  origins  of  the  critical  Reynolds  number  for  the  onset  of  drag  reduction  can  be 
understood  by  analyzing  the  relevant  lengthscales  in  the  flow.  If  the  drag  reduction  and  the  slip  length 
were  dependent  on  the  microridge  geometry  and  channel  dimensions  alone,  as  is  the  case  in  laminar 
flows,  then  we  would  expect  to  find  the  drag  reduction  and  slip  length  to  be  independent  of  Reynolds 
number.  In  turbulent  flows,  however,  there  is  a  third  lengthscale  of  importance,  the  thickness  of  the 
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Figure  9.  Drag  reduction  as  a  function  of  Reynolds  number  for  a  channel  with  a)  a  single 
superhydrophobic  wall  w=30jim  rf=30pm  (•)  and  b)  two  superhydrophobic  walls  containing  w=30|xm 
and  rf=30pm  microridges  (o)  and  >v=60|im  and  c^=60jxm  microridges  (■). 


viscous  sublayer  which  extends  out  toy^  =  5  .  Although  the  viscous  sublayer  thickness  remains  fixed  in 
wall  units,  in  dimensional  form  the  thickness  of  the  viscous  sublayer  decreases  with  increasing  Reynolds 
number  as  y^^^  =  5vyjp/ .  Close  to  the  wall,  where  viscous  stresses  dominate,  the  analytical  solutions 

of  Philip  [31,  32]  show  that  the  influence  of  the  shear-free  air-water  interface  extends  to  a  distance 
roughly  equal  to  the  microridge  spacing,  w,  into  the  flow.  Thus  for  the  superhydrophobic  surface  to 
impact  the  turbulent  flow,  the  microridge  spacing  must  approach  the  thickness  of  the  viscous  sublayer, 
w  «  y^i ,  or  in  other  words  =  3;^  «  5  .  The  microfeature  spacing  in  wall  units  is  at  least  >  0.75  for 

all  the  30-30  surfaces  tested  and  >  2.4  for  the  60-60  surfaces.  The  values  are  calculated  from  shear 
stress  measured  at  the  superhydrophobic  surface.  This  means  that  the  microfeature  spacing  is  minimally 
15%  to  50%  of  viscous  sublayer  thickness  almost  immediately  after  the  turbulent  transition.  Hence  for 
30-30  and  60-60  ridges,  drag  reduction  is  noticed  almost  as  soon  as  a  turbulent  flow  develops.  In  laminar 
flows,  significant  drag  reduction  is  noted  at  feature  to  height  ratios  comparable  to  those  seen  with  the 
present  feature  size  and  viscous  sublayer  thickness  [33].  A  similar  scaling  has  been  observed  for 
turbulent  flow  over  wetted,  rough  surfaces,  where  the  effects  of  roughness  are  not  observed  until  the  size 
of  the  roughness  exceeds  the  thickness  of  the  viscous  sublayer  [34].  As  the  Reynolds  number  increases 
and  the  thickness  of  the  viscous  sublayer  is  further  reduced,  the  presence  of  the  superhydrophobic  surface 
will  more  strongly  influence  the  velocity  profile  within  the  viscous  sublayer  and  reduce  the  momentum 
transferred  from  the  fluid  to  the  wall  and  the  vorticity  of  the  fluid  at  the  edge  of  the  viscous  sublayer. 
Turbulence  intensity  is  thereby  reduced,  increasing  the  drag  reduction.  One  therefore  expects  that 
saturation  of  the  turbulent  drag  reduction  is  likely  in  the  limit  of  very  large  Reynolds  numbers  where  the 
microridges  are  much  larger  than  the  viscous  sublayer.  In  this  limit,  the  drag  reduction  should  approach  a 
limit  of  Df^  =  w/(d  +  w)  as  momentum  is  only  transferred  from  the  solid  fraction  of  the  superhydrophobic 
surface  and  the  viscous  sublayer  is  thin  enough  that  the  no-slip  and  shear-free  portions  of  the  surface  can 
be  considered  independently.  For  the  present  shear  free  area  ratios,  this  limit  would  be  50%.  This  is 
consistent  with  both  the  asymptotic  value  of  our  PIV  and  pressure  drop  measurements.  Drag  reduction 
results  shown  in  Figure  9  appear  consistent  with  this  hypothesis,  the  60-60  ridges  already  appearing  to 
plateau.  As  the  critical  Reynolds  number  will  decrease  with  increasing  feature  spacing,  coarser 
superhydrophobic  surfaces  will  begin  to  perform  better  at  lower  Reynolds  numbers.  It  is  therefore 
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Figure  lOVelocity  profiles  (a)  and  Reynolds  stress  profiles  (b)  predicted  from  DNS  of  channel  flow  at 
a  friction  Reynolds  number  of  Re^  =180  over  superhydrophobic  surfaces  with  microridges  having 

width-spacing  of  15///w-15///w  (□),  30//7w-30//ot(A),  30//7w-50//m(0),  and  30///w-90//ot(V). 

Additionally,  smooth  channel  flow  is  shown  for  reference  (...).  The  symbols  are  employed  to 
differentiate  lines,  and  do  not  reflect  actual  data  points. 

expected  that  equivalent  drag  reduction  performance  will  be  accessible  to  much  finer  microfeature 
spacings  at  higher  Reynolds  numbers.  With  fine  superhydrophobic  surfaces,  little  drag  reduction  may  be 
evident  until  the  viscous  sublayer  shrinks  significantly,  well  past  transition.  This  result  appears 
promising  for  possible  commercial  applications  of  this  technology.  This  is  because  small  feature  spacing 
results  in  a  more  robust  superhydrophobic  surface  capable  of  maintaining  a  coherent  air-water  interface 
at  larger  static  pressures,  while  at  the  same  time  ships  that  might  benefit  from  such  surfaces  operate  at 
Reynolds  numbers  significantly  greater  than  those  interrogated  in  the  present  experiments. 

Numerical  Simulations 

The  DNS  code  has  now  been  extensively  validated  for  turbulent  flows  against  the  benchmark 
simulations  of  Kim,  Moin  and  Moser  [35]  for  grid  resolutions  of  128^  ,  256^  and  512^  having  /?e,  =  180, 
395  and  590  respectively.  The  quantitative  agreement  between  our  code  and  the  benchmark  simulations 
gave  us  the  confidence  to  simulate  ultrahydrophobic  surfaces  having  various  arrangements  of 
microridges  and  microposts  at  a  number  of  different  Reynolds  numbers.  As  one  example  of  the  results 
from  our  DNS  simulations.  Figure  lOontains  the  average  velocity  profile  and  turbulent  Reynolds  stress 
profiles  for  an  ultrahydrophobic  surface  containing  a  series  of  microridge  geometries.  The  simulations 
predict  a  large  slip  velocity  along  the  ultrahydrophobic  wall,  UsUplUmean  >  75%,  and  a  shifting  of  the 
velocity  towards  the  ultrahydrophobic  wall.  Additionally,  the  Reynolds  stresses  show  a  reduction  in  the 
Reynolds  stresses  near  the  ultrahydrophobic  wall,  Ar>  35%,  and  a  non-zero  Rn  and  R33  due  to  the  slip 
velocity  in  both  the  flow  and  vorticity  direction  along  the  ultrahydrophobic  wall.  In  Figure  11,  the  slip 
velocity  and  shear  stress  reduction  is  compiled  for  all  of  the  simulations  performed  to  date  at  Re,  =  1 80. 
The  trends  clearly  show  an  increased  effectiveness  with  increased  microridge  or  micropost  spacing. 
Additionally,  the  simulations  show  that  for  a  fixed  spacing  that  microposts  outperform  microridges, 
however,  if  we  compare  the  two  on  the  basis  of  percentage  of  shear-free  air-water  interface,  then  the 
microridges  are  more  effective.  In  our  most  recent  simulations,  we  found  that  that  unlike  in  laminar 
flows,  in  turbulent  flows,  increasing  the  flow  rate  and  Reynolds  number  results  in  an  increase  in  the 
superhydrophobic  drag  reduction  and  slip  velocity,  although  the  slip  length  appears  to  go  through  a 
maximum  before  decreasing  with  increasing  Reynolds  number  at  high  Reynolds  numbers. 
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Figure  10.  Slip  velocity  and  shear  stress  reduction  for  a  series  of  different  ultrahydrophobic 
surfaces  at  a  Reynolds  number  of  Re^  =  1 80 
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